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Abstract 



We investigate the master nonlinear partial differential equation that governs the evo- 
lution of shear-free spherically symmetric charged fluids. We use an approach which 
has not been considered previously for the underlying equation in shear-free spherically 
symmetric spacetimes. We derive a fifth order purely differential equation that must be 
satisfied for the underlying equation to admit a Lie point symmetry. We then perform 
a comprehensive analysis of this equation utilising the Lie symmetry analysis and direct 
integration. This enables us to reduce the fifth order equation to quadratures. Earlier 
results are shown to be contained in our general treatment. 
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1 Introduction 



The Einstein-Maxwell system of equations plays a crucial role in relativistic astrophysics 
when describing spherically symmetric gravitational fields in static manifolds. In these 
situations we are modeling charged compact objects with strong gravitational fields such 
as dense neutron stars. Several recent treatments, including the works of Ivanov [Ij and 
Sharma et al [2], demonstrate that the presence of the electromagnetic field affects the 
values of redshifts, luminosities and maximum mass of a compact relativistic star. The 
presence of electric charge is a necessary ingredient in the structure and gravitational 
evolution of stars composed of quark matter as illustrated in the treatments of Mak 
and Harko [3] and Komathiraj and Maharaj [4j. Therefore exact models describing the 
formation and evolution of charged stellar objects, within the context of full general 
relativity, are necessary. Lasky and Lun [5l [6], amongst others, studied the role of 
electromagnetic fields in gravitational collapse, formation of black holes and the existence 
of naked singularities. Significant electric fields are also present in phases of intense 
dynamical activity, in collapsing configurations, with time scales of the order of the 
hydrostatic time scale for which the usual stable equilibrium configuration assumptions 
are not reliable (as shown in the treatments of Di Prisco et al [7\ and Herrera et al [8]). 

When seeking exact solutions to the Einstein field equations for an uncharged case, 
the common practice is to make some simplifying assumptions. The assumptions involve 
shear-free and spherically symmetric spacetimes. This approach allows for the reduction 
of the field equations to a single second order partial differential equation. It is re- 
markable that the resulting equation can be treated as an ordinary differential equation. 
For this case the first general class of solutions obtained was a result of the treatment 
by Kustaanheimo and Qvist [9J. Srivastava [10] and Sussman [11] provide extensive 
treatments of the uncharged case. The generalisation to include the electromagnetic 
field is easily performed and is described by the Einstein-Maxwell system of equations. 
Again in this case the field equations can be reduced to a single second order differential 
equation which, interestingly, can also be treated as an ordinary differential equation. 
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Krasinski [12] presents a review of known solutions for the charged case that admit a 
Friedman hmit. A comprehensive study of the physical and mathematical properties 
of the Einstein-Maxwell system in spherical symmetry has been performed by Sussman 
[T3|,[T1] and Srivastava [15]. 

The main objective of this paper is to study the integrability features of the underlying 
differential equation for the Einstein-Maxwell system of field equations. In section 2, we 
reduce the system of Einstein-Maxwell field equations, generalising the transformation 
of Faulkes [16], to a single nonlinear second order partial differential equation. This is a 
governing equation modeling the evolution of shear-free spherically symmetric charged 
fluids. The conventional approach is to use different mathematical techniques to solve 
the master second order equation. The recent treatment of Kweyama et al [T7] is a 
comprehensive attempt to investigate the conditions for the equation to admit a first 
integral or be reduced to quadratures using the Lie symmetry methods of differential 
equations. Here we consider a different approach which we believe is unique and has not 
been considered before for the governing equation in shear-free spherically symmetric 
spacetimes. For the existence of a Lie point symmetry we derive, in section 3, a fifth order 
purely differential equation that must be satisfied. A detailed analysis of this equation is 
performed. We solve, for the first time, the relevant integro-differential equation arising 
from the integration of the fifth order differential equation. A brief discussion of the 
results is given in section 4. 



2 Field equations 

We consider the shear-free motion of a spherically symmetric perfect fluid in the presence 
of the electromagnetic field. We choose a coordinate system = (t, r, 9, 0) which is both 
comoving and isotropic. In this coordinate system the metric can be written as 
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where v and A are the gravitational potentials. We are investigating the general case of a 
self-gravitating fluid in the presence of the electromagnetic field without placing arbitrary 
restrictions on the potentials. For this model the Einstein equations are supplemented 
with Maxwell equations. The Einstein field equations for a charged perfect fluid can be 
written as the system 

P = 3|;-^(2V + A^ + ^)-^ (2a) 
P = J,(-3A?-2A. + 2.A) + Jx(^^ + 2-.A. + ^ + ^)+^ (2b) 
P = ^(-3A?-2A« + 2z.,A,) + ^(i^.. + i^.' + ^ + y + A„) (2c) 

= UrXt-Xtr (2d) 

Maxwell's equations yield 

E = rV-'^^^, = arV^ (3) 

In the above p is the energy density and p is the isotropic pressure which are measured 
relative to the four- velocity = {e^'^, 0, 0, 0). Subscripts refer to partial derivatives with 
respect to that variable. The quantity E = E{r) is an arbitrary constant of integration 
and a is the proper charge density of the fluid. We interpret E as the total charge 
contained within the sphere of radius r centred around the origin of the coordinate 
system. Note that $r = -^lo is the only nonzero component of the electromagnetic field 
tensor Eab = (pb;a — (pa;b whcrc (pa = {^{t, r), 0, 0, 0). The Einstein-Maxwell system 
is a coupled system of equations in the variables p, p, E, a, v and A. 

The system of partial differential equations ([2]) can be simplified to produce an un- 
derlying nonlinear second order equation. Equation (12dp can be written as 

Vr = (InAt)^ 

Then ^ and ^ imply 
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and the potential v has been ehminated. The Einstein field equations (|2]) can therefore 
be written in the equivalent form 

p = 3e''-e-''(2Kr + Xl + ^)-^ (4a) 



= Xte-'' (4c) 
eM A„. - A.^ - ^"l = -F-^ (4d) 



for a charged relativistic fiuid. In the above h = h{t) and F = F{r) are arbitrary 
constants of integration. Equation ( l4dl) is the condition of pressure isotropy generalised 
to include the electric field. To find an exact solution of the field equations, we need to 
specify the functions h, F and E and solve equation (I4dp for A. We can then compute 
the quantities p and p from ( l4al) and (14b I) , and a follows from (j3]). 

It is possible to write ( ]4dl) in a simpler form by eliminating the exponential factor e^. 
We use the transformation, first introduced by Faulkes [16j for neutral fiuids, which has 
the adapted form 

X = r"^ 

y = e-^ 



fix) 



F 

47-2 
2^6 



Then f l4d|) becomes 



y" = f{x)y'' + g{x)y^ (5) 



which is the fundamental equation governing the behaviour of a shear-free charged fiuid. 
In (jS]) the primes denote the derivatives with respect to x. Observe that ([S]) is a nonlinear 
partial differential equation since y = y{t, x). When g = then y" = f{x)y'^ for a neutral 
fiuid which has been studied by Maharaj et al |T8] amongst others. 



6 



3 Lie analysis 



We can verify that 

d d 
G = a{x)— + ib{x)y + cix))— (6) 
ox oy 

is a symmetry of (E]) and the relationship among the functions a{x), b{x), c{x), f{x) and 
g{x) is given by the following system of ordinary differential equations 

a" = 2b' (7a) 

b" = 2/c (7b) 

c" = (7c) 

af' + {2a' + b)f = -3cg (7d) 

ag' + {2a +2b)g = (7e) 

We now combine equations (!7aj) - (!7e|) into one fifth order ordinary differential equation. 
From (!7a|) we have 

2b = a' + a (8) 

where a is an arbitrary constant, from (17b[) we have 

a'" 

and from (17c|l we have 

c = Co + Cix (10) 

Previous analyses (Maharaj et al [18] and Kweyama et al [T7]) have then solved fl7dl) and 
(17e|) for / and (7, respectively. Taking ([9]) into account, Kweyama et al obtained the 
fourth order integro-differential equation for a: 



caa^'"') + 



5a' a > , 

c I 1 — — ca 

'2 2 ' 



a'" = -Ug^c'a-^ exp ( - / ^ ) (11) 



The nature of this particular integro-differential equation has made it difficult to deduce 
much information in general. Here, we focus on obtaining a fifth order purely differential 
equation. 
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We substitute (Ml and (Ml into (iTdl) to obtain 



aa(*^) aa'"c' 5a'a"' aa'" 
We then use ([8]) and (fT2l) in (!7el) to obtain the fifth order differential equation 

2 2 III o 2 I III T r 2 i2 III , A I III , ^ o i I ill r 2 l2 III 

—a c a — Qac a a — 15c a a + Aaacc a + loaca c a — oa c a 
-5ac2a"a"' - Mac^a^'''^ - ISac^a'a^™) + Ga^cc'a^^'') - 2a^c^a^'''^ = (13) 

where c{x) is given by (ITU]) . This is the first time that the symmetry analysis of (^ has 
reduced to solving a fifth order ordinary differential equation. Given that the equation 
is purely differential, one can make recourse to numerical techniques if analytic solutions 
are elusive. 

The fifth order equation (IT^ can be transformed into autonomous form via the 
transformation 

X = A - (lA) 

We end up with the following equation: 

2^2^(^) ^ i3^^'^(») _ 3aAA(™) + 5AA"A"' + 15A'^A"' 

-SaA'A'" + a" A"' = (15) 

Once we obtain solutions to this equation, we can find / and g through direct substitution 
into ^ and ffT^ respectively, after inverting flT^ . 

We thus have the result that the equation 



a'" r, f CLCL 



(*^) aa"'c' 5a' a'" aa"'\ 



^ 4c ^ V 12c2 ^ 12c3 24c2 24c2 J ^ 



admits the Lie point symmetry 

G = ci{x)— + {{a'{x) + a)y/2 + c(x)) — 

where 

c(a;) = Co + Cix 



8 



and a{x) satisfies the fifth order ordinary differential equation 

2^2^(t,) ^ i3^^'^(») _ 3a^^(™) + 5AA"A"' + 15A''^A" 
-SaA'A'" + a" A'" = 

with 

X ^ A ^^'^ 



As observed previously [T7], this is a general result incorporating the ad hoc approach 
of Kweyama et al [19] and Wafo Soh and Mahomed [20]. However, note that the fifth 
order equation f|T5l) has not been derived previously. 

If we let a = in (1151) we obtain 

2^2^('') + 13AA'A(^^) + ^AA"A"' + ISA'^A'" = (16) 

In carrying out the Lie analysis of flT6|) . using PROGRAM LIE, we find that it has the 
following three symmetries 

- w 

We use V\ to reduce the order of (ITBl) . The variables for reduction are 

u = A, V = A 

and the reduced equation is 

2m\S(*") + UuWv'" + Uuv\"' + ^uW^ + 22u\^v%" + bluv^v'v" 

+lbv\" + 2u^vv'^ + lSuv%'^ + 15t;^t;'2 = (17) 

We analyse (|T7|) for symmetries, and obtain the following three symmetries using PRO- 
GRAM LIE: 

ou 



d 

X2 = V — 

ov 



X3 = 2u'^^ + uv 



— + UV — 

du dv 
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The symmetry X3 determines the variables for reduction of ([T 



r = u ^^/^^t;, s = u^^'^v' - -u^^'^v 

2 

and the reduced equation in terms of new variables r and s is 

2r^sh"' + Sr^s^s's" + Ur^s^s" + 2r^ss'^ + 22r^s^s'^ + 22r^sh' + 2rs^ = (18) 
The Lie symmetry analysis of (fT8|) gives three symmetries, namely 

= 4 

or 
d 

OS 

^ d s d 

3 dr r ds 

The reduction variables generated by f/3 are 

p = rs, q = rs' + s 

and the reduced equation is 

2p^q'^q" + 2p^qq'^ + Sp'^q^q' + 2pq^ = 

which simplifies to 

p^qq" + p^q''^ + Apqq' + = (19) 
Note that equation (IT^ is in fact 

{pYY = (20) 

and the solution of fl2Ul) is 

g^ = ^ + ^ (21) 
p 

Note that this result contains the result obtained by Kweyama et al [17]. Usually, when 
an nth order differential equation admits an m-dimensional Lie algebra of symmetries, 
with m < n, there is little hope for the solution of the equation via those symmetries. 
However, in this case we have been able to reduce the equation to quadratures due to 
the existence of hidden symmetries 
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We now return to (fT3|) in the event that a 7^ 0. Integrating the equation ( IT3|) once 
yields the fourth order integro- differential equation 

adx^ 



5a' a 



c'a 



■12g2C^a ^exp 



(22) 



Integrating fl22|) yields 
1 ,„ 



'2a' 



exp 



~2a~ 



dx 



(23) 



where /2 and g2 are arbitrary constants of integration. We multiply f l23|) by a and then 
integrate to get 



aa"'dx = 2/2/ — 6g2 



ca ^^/^^exp 



ca ^^/^^exp 



adx 



2a 



dx 



adx 
~2a 

dx + M 



where M is an arbitrary constant of integration. From (12^ we have 

M = -aa" - la'^ - 2/2/ + 3g2l^ 
2 4 



(24) 



(25) 



where 



ca "^^/^^ exp 



adx 
~2a 



dx 



Again we multiply fl2S]) by al and then integrate to obtain 



aa"'Idx = 2/2 / ca 



~(3/2) 



exp 



6(?2 / ca exp 



adx 
~2a 
adx 
~2a 



Idx 
I^dx - N 



(26) 



where is an arbitrary constant of integration. After performing the integrals in fl26|) 
we find that 



A^ 



-a-(i/2)g^p 
'1 



adx 



It ^12 ^ 



2a 
2 



/ 1 , 1 

ac a c H — ac 

2 2 



-aa --a" + —]I + f2l'-2g2l' 



(27) 
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From fl25l) we have 



^aa" - = M + 2/2/ - 3g2p (28) 



When substituting fl28|) in f l271) we obtain the following equation 

;V . -„-a/.e.p(-/f^)(<.,-i.c.i„e)-(M4)/ 

-f2l' + 92l' (29) 



We thus have the result that the equation 



„ _ a'" 2 f aa'^'"'' aa"'c' 5a' a'" aa"'\ 3 
admits the Lie point symmetry 

d d 
G = a{x)— + {{a'{x) + a)y/2 + c(x))— (30) 
ox oy 

where 

c(x) = Co + Cix 

and a{x) is constrained by the first order integro-differential equation 

N = -a-(^/^) exp (- I ^) (ac' - ^a'c + ^ac^ - (^M + ^) / - ^/^ + g,l' 

and M, A^, /2 and (72 are all arbitrary constants of integration. This is the first time that 
the admittance of symmetry by ([5]) has been reduced to solving essentially a first order 
equation. This is a general result incorporating both the ad hoc approach of Kweyama 
et al [12] and the Noether symmetry results of Wafo Soh and Mahomed 



In the event that a = 0, we have that f l29|) reduces to 
this time with 

' ca-^^/^^dx 



(31) 



We can now integrate fl3ip and obtain the same result as in Wafo Soh and Mahomed 
j. See Kweyama et al p/Tj for a full discussion of this case. 
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It remains to solve (129|) when a ^ 0. If we set 

Ci , aCf 



X 



Co + Cix 



A 



then (129]) reduces to 



N = Ci 



_^-{i/2) 



exp 



+ 



a 



{Co + c,xy 



a 



(32) 



C,{M+^jI-Clhl'-C',g2p 

(33) 



where, in this case, 



^-(3/2) 



exp 



dX 



We can now integrate (1551) further to obtain 



^1/2 



exp 



+ P (34) 



which is an implicit solution for ( IT3|) with a 7^ 0. We note that no reductions of (129|) (or 
its fourth order counterpart) have been previously found for nonzero a. 

Having found a via ( IMl) and (l32l) we now focus on ([5]). Using (l30l) we can transform 
([5]) into the autonomous form 

„2~ 



a 



Y" + aY'+(M+^]Y = + g2Y^ + iV 



(35) 



via the transformation 

X = I — , Y = yexp 
a 



bdx 



exp 



bdx 



dx 



a / J a \ J CL 
The form (!35|) can be more directly analysed. However, when a 7^ 0, we find that we 
cannot directly reduce ( 15B]) to quadratures. Further (Lie) analysis yeilds that, when 
/2 7^ 0, (72 7^ then (|35|) has the following two symmetries 



- dx 



Go 



d_ 
dX 



W^9gJdY 



(36) 
(37) 



provided the following conditions are satisfied 



/I 

M = -— . 

3g2 36 



AT- 2aV2 



27^71 27^?2 



(38) 
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Here, G2 is a new symmetry - Gi is just the transformed form of (l30l) . If we further 
transform using G2 we obtain 

y" = 92y' 

with solution 



where 



dy 



"fy^+c 



a V 35(2/ 

As noted previously [T7] , the values in fl38p correspond directly to a simplification of the 
eigenvalue problem associated with a dynamical systems analysis of (!35|) . 



4 Discussion 



We have reduced the system of Einstein-Maxwell field equations to a single nonlinear 
second order partial differential equation modeling the behaviour of spherically symmet- 
ric charged fluids. This was achieved by utilising the generalised transformation due to 
Faulkes [TB]. The master equation is a partial differential equation because y = y{x,t). 
However it may be treated as an ordinary differential equation. We then derived (for 
the first time) a fifth order nonlinear differential equation. The solutions of this fifth 
order equation yield the solution to the governing Einstein-Maxwell field equations. The 
equation was then transformed into an autonomous form. In the case of a 7^ 0, we 
performed a Lie symmetry analysis of the autonomous fifth order differential equation. 
This allowed us to find a solution to the equation. On letting C2 = 2 in ( |2T1) . we regained 
the result obtained by Kweyama et al [IT] . In our approach we have eliminated the need 
to solve a nonlinear integro-differential equation which was a necessary feature of earlier 
analyses. 

For the a 7^ case, we were not able to make progress via a Lie symmetry analysis 
of the autonomous fifth order equation. However we were able to integrate the non- 
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autonomous fifth order differential equation directly. We obtained an implicit solution 
for this case - a solution that has not been obtained previously. 
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